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ABSTRACT 

We report the discovery of pulsed X-ray emission from the compact source 
IE 1841—045, using data obtained with the Advanced Satellite for Cosmology and 
Astrophysics. The X-ray source is located in the center of the small-diameter supernova 
remnant (SNR) Kes73 and is very likely to be the compact stellar-remnant of the 
supernova which formed Kes73. The X-rays are pulsed with a period of ~ 11.8 s, and a 
sinusoidal modulation of roughly 30%. We interpret this modulation to be the rotation 
period of an embedded neutron star, and as such would be the longest spin period for 
an isolated neutron star to-date. This is especially remarkable since the surrounding 
SNR is very young, at ~ 2000 yr old. We suggest that the observed characteristics 
of this object are best understood within the framework of a neutron star with an 
enormous dipolar magnetic field, R~8xlO^^G. 


Subject headings: pulsars: individual (IE 1841—045) — stars: neutron — magnetic 
fields — (ISM:) supernova remnants: individual (Kes73) 


1. INTRODUCTION 

In 1934, Baade and Zwicky published a prophetic paper making a phenomenological 
connection between supernovae (SNs), the core-collapse of massive stars, and the formation 
of neutron stars (then hypothetical); all purely on grounds of energetics. Decades later their 
conjecture was first vindicated by the discoveries of young pulsars in the Crab and Vela supernova 
remnants (SNRs), and now in a handful of other Galactic SNR. 

Supernova remnants come in at least two distinct morphological types, i.e., shells and plerions 
(Weiler & Sramek 1988), and a majority are the result of core collapse in massive progenitors 
(the non-Type la SNs; van den Bergh &: Tammann 1991). The Baade-Zwicky picture, in its 
simplest interpretation, is somewhat problematic. A majority of SNR appear not to contain either 


^Also: Universities Space Research Association 



- 2 - 


central pulsars or pulsar plerions (as pulsars are beamed, plerions ought to be more commonplace 
than pulsars in the interiors of shells). The predominant “hollowness” of shell-remnants is 
ill-understood, and poses questions about the fate of most core collapses of massive stars. This 
conundrum is nowhere more apparent than in the studies of the youngest SNR, especially those of 
the historical supernovae (Strom 1994). Of the eight historical supernovae which have expanded 
into full-blown SNR, only the Crab Nebula has a pulsar. There is weak plerionic activity (but 
no beamed pulsars) in two others (SN 386 A.D. and SN 1181 A.D.; see Vasisht et al. 1996). It 
follows, therefore, that there is a need to give up our notions about the birth properties of young 
neutron stars, best typified by the Crab pulsar; fast rotation ( .:$ 0.1 s) and a dipole field strength 
clustered around 3 x 10^^ G. That neutron stars may be born in a fashion drastically different 
from the Crab has become increasingly evident via recent X-ray studies. Preliminary evidence of 
this kind includes; the discovery of radio-quiet, cooling neutron stars (Vasisht et al. 1997 and 
refs, therein; Gotthelf, Petre & Hwang 1997) in SNR, the association of the exotic soft gamma-ray 
repeaters with young (~ 10^ yr old) SNR (see Thompson & Duncan 1995), and observations 
of magnetically dominated plerions (Vasisht et al. 1996). Also, the slowly-spinning (P ~ 7 s), 
anomalous X-ray pulsar (AXP) in the 10^ yr old SNR CTB 109, has been known for several 
years (Gregory & Fahlman 1980), although its nature is still widely debated (Mereghetti & Stella 
1995; van Paradijs, Taam & van den Heuvel 1995; Thompson & Duncan 1996). 

This Letter discusses IE 1841—045, an unresolved Einstein point-source discovered near the 
geometrical center of the shell-type SNR Kes73 (Kriss et al. 1985). The rehned Rosat HRI 
location of the object is aj 2 ooo = 18^41™19.2^ and 5j2ooo = —04°56'12.5'' (~ 3" at 90% conhdence; 
(Helfand et al. 1994). The SNR shows no evidence for an extended plerionic core from either radio 
brightness morphology, polarization properties or spectral index distribution. This suggests that 
Kes 73, in spite its inferred youth, lacks a bright radio plerion. To date, no optical counterpart to 
the central X-ray source has been identified for IE 1841—045. Herein, we present the discovery 
of ~ 12 s pulsed X-rays from IE 1841—045 and argue that the source is young and unusual. 

In our companion paper (Gotthelf & Vasisht 1997; hereafter GV97) we present the results of 
imaging-spectroscopy of Kes 73 and the compact source. 


2. OBSERVATIONS 

Kes 73 was observed with the ASGA Observatory (Tanaka, Inoue &: Holt 1994) on 1993, 
October 11-12, as a Performance and Verification (PV) target. Data were acquired by the two gas 
imaging spectrometers (GIS2 and GIS3) and collected with a photon time-of-arrival resolution of 
488 /is in medium bit-rate mode and 64 jxs in high bit-rate mode. We used data made available in 
the ASGA public archive, screened with the standard REVl processing to exclude time intervals 
corresponding to high background contamination, i.e., from Earth-block, bright Earth, and SAA 
passages. An effective exposure of ~ 3.5 x 10^ s was achieved with each detector and the on-source 
measured count rates for the GIS2 and GIS3 instruments were 1.55 (GIS2) and 1.66 (GIS3) counts 
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s“^, respectively. Here, we concentrate on the GIS data exclusively and present data from the two 
Solid-state Imaging Spectrometers (SIS) on-board ASCA, in our companion paper (GV97), which 
reports on spectral and imaging results. We summarize the spectra pertinent to this paper below. 

The spectrum of IE 1841—045 is fit by an absorbed, soft power-law of photon index 
r ~ 3.0 ±0.2 {S,y oc The foreground absorption towards Kes73 is found to be Nh ^ 2 x 10^^ 

cm“^, and is consistent with the kinematic distance estimate of 7 kpc (Sanbonmatsu & Helfand 
1992). The power-law spectral normalization is found to be consistent with no long term spectral 
variation when compared to the count rates observed by Rosat HRI, ~ 0.02 cps (Helfand et al. 
1994). We deduce an unabsorbed model flux of 6.3 x 10“^^ erg cm“^ s“^ (0.5 - 10.0 keV) yielding 
a source luminosity. Lx — 3.5 x 10^^df erg s“^; the SNR distance is T.Odr kpc. 

The long term temporal variability of IE 1841—045 has been examined by Helfand et al. 
(1994) who found no concrete evidence for flux changes on the 10 year baseline between Einstein 
and Rosat. We examined the ASCA data for aperiodic variability by selecting source photons 
from a 5' radius aperture and binning them on ~ 96 min (ASCA orbital period), 10 min and 1 
min durations. The obtained light curves were tested against a uniform model, but no evidence 
for significant variations on these time-scales was found. 

A search for coherent pulsation from the central object was made by combining the two GIS 
high time resolution datasets (tg = 488 /rs). Photons were selected from the entire SNR region of 
~ 5', centered on the compact object from (i) the entire energy band (0.5 - 10.0 keV) (ii) the hard 
band (2.5 - 10.0 keV). The photon time-of-arrivals were barycentered and binned with resolutions 
of 488 fis, 32 ms and 0.5 s. The barycentering and binning procedures were tested on a series of 
datasets of the Crab pulsar and PSR 0540—69. Pourier transforms on the entire datasets were 
performed at each time resolution, interesting periodicities were harmonically summed and later 
folded at the period of interest. 

A significant high-Q X-ray modulation with no overtones is seen at a period P ~ 11.766684 
s (/ = 0.08498571 ± 0.00000004 Hz; see Figs 1 and 2). The modulation is obvious in all the 
above datasets and separately in either GIS on-source time series; it is not observed in off-source 
GIS data, making it unlikely to be an instrumental artifact. The period emerges with greatest 
significance in time-series which contain emission mainly from the central source (hard-band, 

2.5 - 10.0 keV), and is not significant in the soft energy band between 0.5 - 2.5 keV, where the 
nebula is dominant. We suggest that the central object is weakly pulsed at 11.7667 s, possibly a 
neutron star spin period. The pulsed luminosity is Lx ~ 5 x 10^^erg s“^, and the modulation 
level is about 35% of the steady flux from the compact source after subtraction of the estimated 
contribution of the background and the SNR thermal component above 2.5 keV. 

With a period in hand, we have reanalyzed the 18-ks of Rosat HRI data obtained between 
March 16-18, 1992 (Helfand et al. 1994), selecting the few (~ 650) source and background photons 
available from the vicinity of IE 1841—045. We perform a conditional search in a small range 
of periods around 11.76 s using the folding technique (12 phase bins per fold); the resulting 
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periodogram is displayed in Fig lb. We cautiously forward the suggestion of a peak-up at a 
barycentered period of 11.7645 s (4.0 cr; which is the expected significance given the HRI count rate 
and the ASCA derived modulation); the observed FWHM of the periodogram excess of 10“^ s 
is consistent with expectations, given the period and the time-span of the HRI observations. Also, 
the crest and trough in the HRI profile (Fig 2b) match those of the GIS profile quite well. A linear 
interpolation, assuming steady spin-down, gives a period derivative of P ~ 4.73 x 10“^^ s s“^. 


3. DISCUSSION 

Our interpretation of IE 1841—045 is that of a young neutron star that was born during a 
supernova that now forms the SNR Kes73. The kinematic distance to the SNR of 6.7 - 7.5 kpc is 
consistent with its high foreground X-ray absorption, Nh — 2 x 10^^ cm“^. The small shell radius 
{Rg — 4.7^7 pc), along with intense radio and X-ray shell emission are characteristics of a young 
supernova remnant. This notion is supported by our detection of a hot thermal {kT ~ 0.8 keV) 
X-ray continuum in the shocked gas. Helfand et al. (1994) have argued that the SNR is likely 
to be in transition between free expansion and the adiabatic phases, whereby, the Sedov age of 
Tg < 2500 yr must be an upper bound. The SIS spectra show the enhancement of Mg over the 
species S, Si and Fe, along with possible evidence for highly absorbed emission from O and Ne 
in the 0.5-0.9 keV range (see GV97). During the course of their evolution, massive stars produce 
large quantities of 0-group (O, Ne and Mg) which are ejected during the supernova (see Hughes 
et al. 1994). Hence, there is spectral evidence that Kes73 is young and still ejecta dominated. On 
spectral grounds we also favor a Type H or Ib origin for Kes73, i.e., from a massive progenitor. 
We consider it unlikely that the neutron star was born in an accretion-induced collapse of a heavy 
white dwarf (Lipunov & Postnov 1985). 

The pulsar in IE 1841—045 has common properties with the peculiar X-ray pulsar, IE 
2259-1-586 (Gregory & Eahlman 1980; Gorbet et al. 1995), and the soft gamma-ray repeaters 
(Thompson & Duncan 1996). IE 2259-1-586 is a 7-s spin rotator, and coincides with a 10^ 
yr-old SNR (CTB 109). Much like IE 1841—045, it has a soft X-ray spectrum best represented by 
a blackbody at 0.45 keV with a non-thermal tail with Si, oc and — 0.5 — lx 10^^ erg s“^. 
It has a history of nearly steady spin-down and no detected binary modulation (Iwasawa, Koyama 
& Halpern 1992), optical companion, or quiescent radio emission. There are four other X-ray 
pulsars, 4U0142-I-614 (Hellier 1994; Israel et al. 1994), IE 1048.1-5937 (Seward, Gharles & Smale 
1986), RXJ1838—0301 (Schwentker 1994; possibly also associated with a 10^ yr old SNR), and 
IRXS J170849.0—400910 (Sugizaki et al. 1997) which have low luminosities {Lx ~ 10^® — 10^® erg 
s“^), periods of order 10-s that are steadily increasing, soft spectra and no detected companions or 
accretion disks. In all the above cases, spectra can be fit with soft power-laws with indices in the 
range 2.3 - 3.5 (Corbet et al. 1995). Collectively, these have been grouped into a class called the 
breaking X-ray pulsars (Mereghetti & Stella 1995) or alternatively, the anomalous X-ray pulsars 
(van Paradijs et al. 1995). We observe AXPs through a substantial distance in the Galactic disk 
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(with foreground columns densities ~ 10^^ cm ^), which suggests that they are not commonplace. 

The rotational energy of IE 1841—045 is far too small to power its total X-ray emission of 
Lx ~ 4 X 10^^^ 7 - The maximum luminosity derivable just from spin-down is 

where I is the moment of inertial of the neutron star, and Tg 2 X 10^ yr is the SNR age. The 
mechanisms for powering the X-rays could then be either (i) accretion from a high-mass X-ray 
binary (Helfand et al. 1994), a low mass companion (Mereghetti & Stella 1995), or a fossil disk 
(van Paradijs et al. 1995), or a merged white dwarf (Paczyncski 1990) (ii) intrinsic energy loss, 
such as initial cooling or the decay of magnetic fields in a magnetic neutron star (Thompson & 
Duncan 1996). 


The strongest argument for accretion as the source of energy is that the inferred accretion rate 
is just that required if the NS is close to its equilibrium spin period Peg, with a field B r\j 1012 G 
typical for young pulsars (see Bhattacharya & van den Heuvel 1991) 


eq 


'5 X 1011 




However, only a pathological evolution scenario involving accretion could bring an energetic dipole 
rotator to its present rotation rate within ~ 2 x 10^ yr. The strongest support for IE 1841—045 
as an accretor will be from future identification of an infra-red (large foreground Ay — 10 mag.) 
companion or an accretion disk. As in the case of other AXPs an infrared counterpart may not 
be easily identified (Coe, Jones & Lehto 1994 and refs, therein). The pulsar has properties that 
may already preclude accretion as a power source: (i) high-mass neutron star binaries with a NS 
accreting from the companion wind sometimes go into low luminosity states with Lx 1035 

erg 

s“^, and have periods in the range 0.07 - 900 s. However, in general they display hard spectra 
(0.8 < T < 1.5) and strong aperiodic variability on all time-scales (Nagase 1989). If the X-ray 
source is indeed a high mass X-ray binary then the inferred P could be the result of orbital 
Doppler effects, (ii) Neutron stars with disk accretion from a low mass companion or a fossil 
disk, with the latter having formed from SN debris or a Thorne-Zytkow phase (van Paradijs et 
al. 1995), should display similar accretion noise in the light-curve, (iii) Other AXPs show near 
steady spin-down on time-scales 10^ — 10^ yr, although this is controversial (Baykal & Swank 1996; 
Corbet & Mihara 1997). The long term torque behavior of IE 1841—045 will only be evident with 
future observations, iv) Einally, accretion models would have to be stretched in order to explain 
the slow rotation period (inside a young SNR) and its associated spin-down time-scale, ^ 2500 
yr. First, it difficult for accretion torques to spin-down a pulsar to 12-s in ~ 10^ yr from initial 
periods Pi ^ 10^ ms unless, of course, the pulsar were born a very slow rotator, which is quite 
interesting in its own right. Secondly, if the pulsar were rotating near its equilibrium period, as in 
the Ghosh and Lamb (1979) scenario, the spin-down time of P/2P ~ 3900 yr is inconsistent with 
the implied accretion rate, M ~ 10“^^ Mq yr“^ (assuming the pulsar has a normal dipolar field 
~ 10^2 c); these usually lie in range 10^ — 10^ yr. 



- 6 - 


A dipolar magnetic field vs. period scaling can also be obtained under the assumption that 
the neutron star is isolated, and has undergone conventional pulsar spin-down from torques 
provided by a relativistic wind, as in the Crab. For dipolar secular spin-down, the implied B^ipoie 
is enormous. 


P = 10 s 


\l/2f ^dipole w 


R 


)\-- 


M 


\-l/2 


3xl03yr' AQis G''15 km' A.4 M© 


Such highly magnetized neutron stars (with dipolar field strengths B ~ 10^^ — 10^® G) or 
“magnetars” have been postulated by Thompson & Duncan (1995 and refs, therein) to explain 
the action of soft gamma-ray repeaters. Magnetars have magnetic flux densities that are a factor 
10^ larger than the typical B ~ 10^^ G fields supported by radio or X-ray pulsars and perhaps 
represent a tail of P-field distribution in NS. After birth, they spin-down too rapidly to be easily 
detectable as radio pulsars, assuming that they are capable of radio pulsar action at all. The 
dipole energy in the star’s exterior, a small fraction of the total magnetic energy, exceeds the 
rotational energy of the NS after roughly 2B^^ yr, where B = lO^^Pis G- Magnetism then quickly 
becomes the dominant source of free energy in an isolated magnetar. 


The derived spin-down age of the pulsar ~ 3900 yr, is consistent with the SNR age 
(ages inferred from the estimator P/2P are larger than the true age as they measure linear 
spin-down).The equivalent dipolar field is Bdipoie — 3.2 x 10^® (PP)^/^ 8 x 10^^ G. There is an 

intriguing possibility that the the pulsar in Kes73 was born as a magnetar ~ 2 x 10^ yr ago, and 
has since spun-down to the long period of 11.7-s due to rapid dipole radiation losses. It could 
be unobservable as a radio pulsar due to period dependence of beaming (Kulkarni 1992); it is 
also possible the radio pulsar mechanism may operate differently or not at all above the quantum 
critical field, Per — 4.4 x 10^^ G. In a magnetar, the observed X-ray luminosity would be driven 
by the decay of the stellar B-field via diffusive processes, which set in at an age 10^ - 10^ 
yr (Thompson &: Duncan 1996; we assume that diffusion of field lines through the crust by 
Hall drift, and the core by ambipolar diffusion occur on time-scales Tuaii — 5 x yr and 

Tamh — 3 X lO^Pj^^ yr, respectively; Goldreich & Reisenegger 1992). Magnetic field decay powers 
the star on average, at a steady rate of (l/6td)P^PjjpQ;g 1036 

erg s ^ for the first decay time, 

td ~ 10^ yr. Decay in the core is likely to keep the stellar surface hot via release of internal 
magnetic free energy, while crustal decay is likely to set up a steady spectrum of Alfven waves in 
the magnetosphere which can accelerate particles to produce the soft non-thermal tail observed 
in the pulsar spectrum (GV97). Such a neutron star may in time (~ 10^ yr) display the soft 
gamma-ray repeater phenomenon (Thompson & Duncan 1996). 

In conclusion, we reiterate that IE 1841-045 is in our estimation a young (~ 2000 yr-old) 
neutron star spinning at an anomalously slow rate of ~ 11.8-s, possibly with very strong torques 
on its rotation. The claim of rapid spin-down is based on a weak detection of periodicity in the 
archival Rosat HRI data, and needs to be urgently tested in future observations. Whatever the 
hnal consensus on IE 1841—045, it is a most unusual and exciting object, the understanding of 
whose nature should make us rethink important aspects about the birth process of neutron stars 
as a whole. 
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Fig. — 1 (Top) A power spectrum of photons from both GIS cameras displayed in the range 11.75 - 
11.78 s. Photons were selected from the hard-band (2.5 - 10.0 keV). The main peak near 11.766684 
s is the putative pulsar period. The powerful side-lobe peaks are separated from the main peak 
by 0.00017 Hz, the ASCA orbital period. (Bottom) A periodogram of vs. period of the HRI 
data, peaks up at roughly 4a as is expected from the GIS profile (even though the energy range is 
different). The peak-up period is 11.7645 s. The total number of HRI counts for this was ~ 800. 
The search was done with twelve bins across the folding period. 

Fig. — 2 (Top) A normalized folded profile of the GIS data (including background), with 12 
bins of resolution, and a folding period of 11.7667 s. The profile is roughly sinusoidal, with about 
~ 35% modulation (after accounting for the background). The start epoch of folding is MJD 
4.9271609362 x 10^. (Bottom) A normalized folded profile of Rosat HRI data, with 12 bins of 
resolution, at a folding period of 11.7645 s. The start epoch of folding is MJD 4.8697189 x 10*^. 







































